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ABSTRACT 

We report on a series of three-dimensional magnetohydrodynamic simulations of active galactic 
nucleus (AGN) jet propagation in realistic models of magnetized galaxy clusters. We are primarily 
interested in the details of energy transfer between jets and the intracluster medium (ICM) to help 
clarify what role such flows could have in the reheating of cluster cores. Our simulated jets feature a 
range of intermittency behaviors, including intermittent jets that periodically switch on and off and 
one model jet that shuts down completely, naturally creating a relic plume. The ICM into which these 
jets propagate incorporates tangled magnetic field geometries and density substructure designed to 
mimic some likely features of real galaxy clusters. We find that our jets are characteristically at least 
60% efficient at transferring thermal energy to the ICM. Irreversible heat energy is not uniformly 
distributed, however, instead residing preferentially in regions very near the jet/cocoon boundaries. 
While intermittency affects the details of how, when, and where this energy is deposited, all of our 
models generically fail to heat the cluster cores uniformly. Both the detailed density structure and 
nominally weak magnetic fields in the ICM play interesting roles in perturbing the flows, particularly 
when the jets are non-steady. Still, this perturbation is never sufficient to isotropize the jet energy 
deposition, suggesting that some other ingredient is required for AGN jets to successfully reheat cluster 
cores. 

Subject headings: galaxies: jets - galaxies: clusters: general - methods: numerical - MHD 



1. INTRODUCTION 

Over the past decade. X-ray observations of galaxy 
clusters have clearly illustrated the immense energies as- 
sociated with interactions between active galactic nuclei 
(AGNs) and their environments. For example, discover- 
ies of X-ray cavities associated with A GN radio -emitting 



'relic bubbles' (e.g . McNamara et a/.|2Q QQ; Fa bian et al 
Fujita a/. l'2UU2rNuls en 'a/ll2QQ2HWise et al. 
2QQ7[ |Sanders et al. 2009 ) demonstrate that radio lobes 
associated with current or previous AGN activity often 
completely displace the surrounding intracluster medium 
(ICM). Likewise, X-ray observ ations of 'ghost cavities' 



history of realistic jet lifetimes and duty cycles. Some at- 
tempts to esti mate AGN intermit tency are based largely 
on energetics. Soker et al (2001), for example, estimate 
jet lifetimes of 10*^ — 10^ years with duty cycles of 5% 
or less to stifle cooling while matching optical and ra- 
dio observations of cooling c ores. Observations of bub- 
bles in the absence of AGN (Mazzotta et a/.| |2002 ) have 



devoid o f radio emissi on {e.g. McNamara et al. 2001 



Mazzotta et al. 



also been used to estimate lifetimes of 10'^ — 10^ years 
and duty cycles of at most 10%. Other approaches focus 
more on timescales associated with observed featur es in 
individual so urces, such a s those conducted by Fabian et| 
oTI ( |2003|), [Forman et o TI (|2005[ ), [jamrozy et a/.| ( |2007p 



and Blanton et al. 



20021) have shown that these cavities one event per 10^ 



pDDD) 
10^ 



ot whom tind a rate of 
Observations of the so- 



maintain some structural integrity even after the radio 
population has aged beyond the limits of detectability. 
The energies associated with these plasma bubbles are 
found to be upward of 10^^ - 10^^ e rg (Birzan et a/., 2004] 
Dunn et a/.|2005[jWise et a/.|20Q7| , suggesting that they 



years. 

called double-double radio sources fSchoenmakers e t al. 
|200 0|), in which two distinct episodes of jet activity are 
clearly visible, have im plied rather sh ort AGN lifetimes 
of ~ 10^ years at most. Groom et al. ( 2004) take a com- 
pletely different approach, using SO clustering to esti- 
mate AGN intermittency, and their result sugge sts AGN 
lifetimes of 10^ - 10^ years. Most recently, Birzan et al. 



[Kaastra et a/.,200f[ " Tamura et a/. 2001 ). This idea is fur- 



ther supported by simulations ot cluster formation and 
ev olution with AGN fee dbac k, such as those conducte d 
by jBruggen et a/.| ( |2005[ ) and [Sijacki fc Springel| (|2006|). 

Observations of detached relic bubbles and evolved 
ghost cavities also suggest that a complete picture of 
AGN interactions with the ICM should incorporate a 
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could play a signiticant role m stifling cooling flows and 

maintaining t he ^ 2 keV temp erature flo or observed in ^ ^ , . 

clusters (e.g. "Pet erson et a/.r 2001; Fabia n et a/.||2001[ (|2009) use a complete sample of clusters to estimate duty 

cycles of at least 60% for A GN outbursts in cooling flow 
clusters. Diehl et al. (2008) also examine a large sample 
of clusters and infer that the observed size distribution of 
X-ray cavities is most consistent with either continuous 
hydrodynamic inflation or discrete production of bubbles 
by magnetically dominated jets. 

Although there is much observational evidence to sup- 
port models of intermittent jets that feed energetic 
plasma bubbles to the ICM, much still needs to be un- 
derstood about the precise nature of energy flow in these 
systems. Simple analytic n iodels , such as those devel- 
oped by Cioffl & Blondin (1992), Kaiser & Alexander 
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( |T997| ), and 
tools for ana! 



Komissarov & Falle (1998), are powerful 
yzing the basic dynamics and energetics in 



these systems, but are not designed to address the details 
of energy transfer to the ICM. Furthermore, these mod- 
els tend to feature arguments based upon self-similarity 
or constant jet luminosities, and so may not apply to re- 
alistic ICM atmospheres or intermittent jets. Likewise, 
numerical studies that f ocused on jet dyna mics in both 
the non-relativistic (e.g. Norman et al. 1982; Hardee 
Clarke"1992» 'Norman"l l9961 |Carvallio 7n?Dea||2Qd^ 
Krause 2003, 2005) a nd relativis t ic regi mes [e.g. Dun- 
can & Hughes 1994; Marti et aL|p 97; Alo y et al.|[T^^ 



Mizuno et al. 2007; Keppens et al . 2008; Mignone et al.' 
2009p made great strides toward understanding how jets 
propagate, but were not necessarily interested in the de- 
tails of energy transfer to the ICM. [Rosen et al.| (|1999|) 
provide a useful analysis of the primary^itferences be- 
tween relativist ic and non-relativistic simulations, point- 
ing out that jets of equal thrust have similar ages and 
efficiencies in both cases, although relativistic jets typi- 
cally develop smaller cocoons than their non-relativistic 
counterparts. 

Recent numerical studies of jet energetics, on the other 
hand, have contributed greatly to the underst anding of 
ene rgy flo w from AGN to their environments. Zanni et 
al. (2005|), for example, conducted an ensemble of two- 
Himensional simulations of jets in detailed cluster envi- 
ronments, finding that a significant fraction of jet energy 
is deposited irreversibly into the ICM. Likewise, our pre- 
vious study o f three-dimensiona l magnetohydrodynamic 
(MHD) jets ( O'Neill et a/.|2005D found that energy trans- 
fer to the ICM was etticient tor both uniform and strati- 
fied environments. Continu ing work in three dimensions. 



Vernaleo & Reynolds ( 2006 ) pointed out that simple hy- 
dro dynamic models of jet feedback characteristically de- 
velop a relatively lossless low-density channel and sub- 
sequently fail to dissipate and distribute e nergy widely 
yh to prevent catastrophic cooling. |Heinz et al' 



sim- 



enough to prevent catastrophic cooling. [Memzj 
([2006), however, attempted to address this issue by 
ulating wobbling hydrodynamic jets in an evolving clus- 
ter atmosphere that served to isotropize the deposited 
energy. Whether the jet and jet-blown cocoon remain 
characteristically stable or break apart or mix with the 



ambient medium (as illustrated recently in Gaibler et al 
[2009, for example) clearly influences the etticiency with 
which jets deliver energy to cluster cores. 

Some recent numerical studies of relic bubble propaga- 
tion also have had implications for energy trans fer and 
transport i n the ICM. Simulations conducted by IRobin- 



son et 



dimensions and 



(2004 ) and Jones fc DeY oung (2005 



Kuszkowski et al. 



(|2007j),|C 



in two- 



O'Neill et al. 



^2009), and Dong !Stone| (|2009p nTtliree dimensions, 
for example, illustrated that rising MHD bubbles lift a 
great deal of ICM material in their wake. This material 
has the potential to mix and dis tribute energy within 
the ICM. A recent paper by jBriiggen et al. (2009) points 
out that this mixing at the bubble interface is further 
enhanced in the hydrodynamic case when turbulence is 
treated with a subgrid approach, which explicitly mod- 
els the energetics and sizes of turbulent eddies that are 
not properly resolved by the numerical grid. This is in 
contrast to the usual method of treating turbulence in 
which the computational grid resolution sets the effec- 
tive scale of dissipation in the system. It is currently 



unclear how t his approach will affect th e case of magne- 
tized bubbles. Heinz & Churazov (2005), taking a differ- 
ent approach, showed that hydrodynamic bubbles could 
explicitly transform wave energy into heat that could be 
dissipated. 

Several groups have also explored the effects of jet 
intermittency on flow dynamics and energetics. The 
earliest ex ample of a simulated r estarting jet was con- 
ducted by Clarke & Burns (1991), who found that sub- 
sequent jets moved quickly through previous jet material, 
forming a new bow shock in the pro cess. U sing a two- 
dimensional hydrodynamics code, [Reynolds et al 
found that plumes from dead radio galaxies lift 
stantial amounts of ICM material as they rose and that a 
large fraction of the total jet ene rgy went into heating th e 
ICM. This was followed by Basson fc Al exa nder[ ([20031) , 
who demonstrated in three dimensions that a buoyant 
jet-blown bubble continues to dredge up material from 
the cluster core for muc h longer than the j et lifetime. 
Also in three dimensions, Omma et al. (2004) found that 
most of the jet energy in an outburst ends up too far away 
from the cluster core. In a subsequent simulation that 
included cooling, Omma fc Binney (2004) found that the 
amount of dense, cool gas immediately above the central 
AGN had a pronounced effect on where most of the jet 
energy was dissipated. While the p reviously mentioned 
work by Vernaleo fc Reynolds ( 2006 ) showed that this ef- 
fect was never sufficient to completely offset cooli ng in a 
simple ICM atmosphere, recent simulations by Brii ggen[ 
fc Scannapieco (2009) suggest that the modeling of sub- 
grid turbulence may alleviate some of these problems. 

Our study combines several of these different investiga- 
tions, exploring the interactions of intermittent jets with 
the ICM. We seek to determine how intermittency affects 
energy flow, asking in what form, where, and how rapidly 
energy is transferred from jet-produced structures to the 
ICM. These simulations specifically examine whether a 
combination of jet intermittency, magnetic stresses, and 
ambient density fluctuations alone can lead to broader 
distribution and dissipation of jet energy, or whether 
something like the more chaotic jet model of Hein z et[ 



al. 



.( [2Q06j , or a very hi gh Reynolds number turbulenc e 
such as that proposed by Briiggen fc Scannapieco ( [2009[ ), 
is required. Furthermore, we seek to understand the ba- 
sic dynamics and morphology of intermittent flows. We 
test how flow structures produced by intermittent and 
extinct jets evolve, how they transport and transfer en- 
ergy, and how they interact with and potentially enhance 
ambient magnetic fields. 

To address these questions, we employ a set of MHD 
simulations of jets in realistic cluster environments. The 
jets are oppositely directed, launching into an ICM with 
full three-dimensional gravity, tangled magnetic fields, 
and density fluctuations designed to mimic cluster sub- 
structure. Section[2]describes the details of the numerical 
methods and the models employed. Section |3] contains a 
discussion of the results and astrophysical implications, 
while Section m lists the main conclusions of this work. 

2. CALCULATION DETAILS 

Our simulations employ a second-order Eulerian total 
variation diminishi ng (TVD) non - relativ istic ideal MHD 



code, described in Ryu fc Jones (1995) and Ryu et al. 
(1998). The code explicitly enforces the divergence- free 
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TABLE 1 

Summary of Simulations of Double Jets in 
Galaxy Clusters 



ID 1 


Jet Intermittency 


Final Age 


ST 


Steady inflow for duration 


^ 59 Myr 


113 


Switch on/off every 13 Myr 


?^ 173 Myr 


126 


Switch on/off every 26 Myr 


^ 104 Myr 


RE 


On for 26 Myr, then off 


^ 144 Myr 


AM 


No jets included 


^ 183 Myr 



All models feature identical computational grids 
(600 kpc X 480 kpc x 480 kpc) and ICM structures. 
ICM model features core density po = 8.33 X 10~^^ 
g cm~^ and pressure Pq = 4.0 X 10~^° dyne cm~^, 
corresponding to a core sound speed cq = 895 km s~^ 
and temperature 2.5 keV. Ambient magnetic fields are 
tangled with a characteristic core value Bq ~ 10 /iG. 

In fully active states, jets feature Mach number 
Mj = 30, corresponding to velocity vj — 0.0895c. Jet 
densities are calculated from fixed pj = r] pQ, with 
?7 = 0.01. Jet magnetic fields are completely toroidal 
with maximum strength Bj = 10 /iG. 

condition for magnetic fi elds through a con strained trans- 
port scheme detailed in Ryu et al. (1998). The numer- 



ical method conserves mass, momentum, and energy to 
machine accuracy. Gravitational energy is handled in a 
fully conservative fashion by including this component in 
the total energy of the gas. Additionally, energy fluxes 
across the boundaries and changes in energy on the grid 
are tracked and output for each energy type to facilitate 
analysis of energy flow (see Sect ion 3.1). Model param- 
eters (discussed in Section 2.2) are selected such that 
the cluster cooling times are longer than the simulation 
times (tcooi > 250 Myr in the cluster cores), allowing us 
to neglect cooling in our simulations. 

All simulations described here are performed on a 
three-dimensional Cartesian grid of full physical dimen- 
sions X = 600 kpc, y = z = 480 kpc. Each zone repre- 
sents one cubic kiloparsec of volume with Ax = Ay = 
Az = 1 kpc. Given that the total jet diameter is 14 
kpc (see Section 2.1), this is sufficient to resolve shocks 
and multiple scales of turbulence within the jet and co- 
coon. Wit h th e exception of one model (AM, described 
in Section 2.1), these simulations feature sets of oppo- 
sitely directed jets emanating from a cylindrical internal 
boundary located in the central region of the computa- 
tional grid {x = y = z = 0). The jet region is updated 
prior to each TVD step, and the settings in this region 
are used to control jet luminosity and intermittency. 

Since the jets are launched from the center of the com- 
putational grid, the outer grid boundary conditions have 
little influence on the evolution of jet/cocoon structures. 
Still, the outer boundaries are crucial in maintaining 
the dynamical stability of the constantly evolving am- 
bient medium. In the outer grid boundaries, we use a 
set of modified continuous boundaries designed to main- 
tain approximate hydrostatic equilibrium in the ambient 
medium while minimizing the impact of small-amplitude 
waves incident upon the boundaries. Specifically, we ap- 
ply the constraints that the sound speed be constant and 
that hydrostatic equilibrium apply to set the density and 
pressure conditions across a given boundary. We find 
that these boundaries maintain the total mass and en- 
ergy of the grid to an accuracy greater than 99.5% in the 



absence of a jet over three grid sound-crossing times, a 
time much longer than any of our jet simulations. 

A passive color tracer Cj, representing the mass frac- 
tion of jet material, is introduced through jet flow to 
identify material that has passed through the jet ori- 
fice. Additionally, we include a population of passive, 
nonthermal, relativistic cosmic ray particles in our mod- 
els. These cosmic rays are transported, injected, acceler- 
ated, and aged in a self-consistent fashion to enable the 
construction of realistic synthetic observations (see, e.g. 



Jones et a/.|[T9 99, Tr egillis a/.||2001a[ [Tregillis etH 



Jones et at.\\W )jy, Ir efflllis et a 
2004, Mendygr al et a/.|in prep ~ 



Here, we describe live simulations in which the initial 
ambient conditions are precisely the same, but the time 
histories of jet activity vary. In the following two subsec- 
tions, we describe the jets and ambient media in detail. 
The physical parameters of each simulation are described 
in Table 1. 

2.1. Bi- directed Jets 

Jet inflow is introduced entirely along the ±x-axis 
without any initial transverse component. The previ- 
ously mentioned internal jet boundary consists of a cylin- 
drical region originating from x = and extending 6 kpc 
along ±x. The resulting jets have uniform cores of radius 
rj = 3 kpc, surrounded by a concentric transition annu- 
lus that smoothly connects to the ambient conditions at 
a radius of 7 kpc. The nominal core density contrast is 
T] = pj/pa = 0.01, and the jets initially enter in pressure 
equilibrium with their surroundings. The jet core mag- 
netic field in this set of models is completely toroidal, of 
the form 5^ = Bo{r/rj) within the jet core, but that 
decreases quadratically to zero outside of the core re- 
gion. This jet field was chosen to be strictly toroidal so 
that it could be initially divorced from the ambient field. 
The fiducial ratio of thermal pressure to magnetic pres- 
sure in the jet (measured at the boundary of the core) 
is 13 ^ 100, and the associated (maximum) inflowing jet 
field strength is Bq = 10 pG. 

For this set of models, the jet luminosity can be ex- 
pressed in terms of the steady jet values of the density 
Pj , velocity Vj , and pressure pj as 



Lj ^ 2.12x10"^^ 



1.90 



erg s 

(1) 

ignoring the relatively minor contributions of magnetic 
and gravitational energy to the jet luminosity. Since it is 
the velocity (and not the luminosity) that goes smoothly 
to zero in the transition annulus, the size of the annu- 
lus enters into the kinetic and thermal energy flux terms 
differently, resulting in the prefactor of 1.9 (see the Ap- 
pendix for details) . This analytic expression is more than 
99% accurate over the lifetime of a steady jet. The first 
jet model features two steady Mach 3 jets (labeled 'ST' 
for 'steady') that correspond to a physical flow speed 
Vj ~ 0.1c and jet luminosity Lj = 6.1 x 10"^^ erg s~^ 
per jet. The other runs in this series of simulations fea- 
ture jets with variable activity. In all cases, the maxi- 
mum injected jet speed and Mach number match those 
of the ST model, and it is worth noting explicitly that 
these assumed jet luminosities and speeds (correspond- 
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Fig. 1. — Volume rendering of flow speed for the ST (upper left), 126 (upper right), 113 (lower left), and RE (lower right) models, after 
the jet disturbances have nearly reached the computational grid boundary. The contrast in the RE image has been enhanced to increase 
visibility. Animations of these quantities as seen from several different angles will be available through the electronic version of this paper. 



ing to a Lorentz factor of 1.004) are consistent with a 
non-relativistic treatment. In two runs, the jets intermit- 
tently switch on or off at regular intervals. In these two 
cases, the jets switch states every 26 Myr (labeled '126') 
and 13 Myr (labeled '113'), respectively. This switch 
is accomplished by ramping the density, pressure, and 
momentum density exponentially over a 3.27 Myr (1.64 
Myr) timescale for the 126 (113) run. The exponential 
ramp changes the density and pressure to a volume- 
averaged value sampled from a spherical region imme- 
diately surrounding the jet, allowing for a smooth tran- 
sition from active jets to a quiescent state. In the fourth 
model (labeled 'RE' for relic), designed to mimic radio 
relic sources, the jets are completely shut down after 26 
Myr and are not restarted. Additionally, there is a fifth 



model (labeled 'AM' for 'ambient') in which the jets are 
never activated. This final model allows us to separate 
the evolution of the ambient medium from jet-driven be- 
haviors to correctly assess the role of the environment in 
these simulations. 

Figure [l] shows volume renderings of flow speed at late 
times for each of the four jet-driven models. In each 
frame, the jets originate near the center of the image at 
an angle with respect to the viewer, with the approach- 
ing jet pointing toward the lower right and the receding 
toward the upper left. All four of the jet-driven simula- 
tions run until wave or shock disturbances from the jets 
reach a grid boundary. This timescale varies from model 
to model depending on how the flow is driven, and the 
total simulation times are given in Table 1. 
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Fig. 2. — ICM density, pressure, and temperature profiles. Values 
of ^0, and Tq are given in Table 1. The labels rd and rc2 
represent the locations of the ICM density core radii. 

2.2. Galaxy Cluster Environments 

In all five models, we employ a single type of ambient 
medium designed to incorporate the essential properties 
of galaxy cluster environments. The gas den sity bas e 
profile is a double- /3 model {e.g., pta fc Mitsudaj[2QQ4 ) , 



similar to those typically used to observationally describe 
clusters, of the following form (see Figure |2|: 



Pa{r) = po 



fl 



/2 



(l + (^)^) 



(2) 



where po = 8-33 x 10~^ g cm"'^ is the density at r = 0. 
This profile features two density "cores" weighted differ- 
ently: Tci = 50 kpc with a weight /i = 0.9 and rc2 = 200 
kpc with a weight /2 = 0.1. These values were chosen 
to be consistent with the best-fit values found through a 
statistical analysis conducted by Ota & Mitsuda (2004) 
on 79 X-ray clusters. For these models, fj = 0.7, which 
generates an asymptotic profile only slightly steeper than 
when r > > rc2 > Tca • On top of this profile, we have 
superposed a Kolmogorov (1941) spectrum of density 
fluctuations with a maximium amplitude of ±0.10pa(^) 
locally. This amplitude of adiabatic fluctuations in den- 
sity corresponds to the range of pressure fluctuations ob- 
served by Schuec ker et al.\ ([2004j) in the Coma cluster. 
This initial spectrum of fluctuations is spatially tapered 
by a Gaussian envelope near the edges of the compu- 
tational grid to impede the development of unphysical 
flows at the grid edges. 

For these simulations, gravity is derived from an NFW 
profile ( Navarro et a/. 1997) , intended to model the un- 
derlying cluster dark matter distribution: 



Ps 



(^)(l + ^)2' 



(3) 




Fig. 3. — Initial tangled ICM magnetic field lines. Arrows on 
field lines indicate direction while intensity refiects field strength 
(light: stronger, dark: weaker). Jets emerge from the region near 
the center of the image. 



which generates a gravitational acceleration: 



9{r) 



^^Gr\^Ps 



In 1 



r 

^drr 



(4) 



where rdm = 400 kpc is taken to be the characteristic 
dark matter core radius. Another free parameter is the 
dark matter density, ps. Rather than deriving this from 
estimated cluster masses, we instead selected the value 
that would produce a reasonable temperature distribu- 
tion (see below). In this case, the temperature profile 
selected corresponds to ps ^ 4.3 x 10~^^ g cm~^. For a 
cluster with a virial radius = 2 Mpc, this leads to a 
virial mass My = 5 x lO^^M©, which is within a factor 



of a few of the mass of the Perseus cluster {e.g. Ettori et 
al. 1998). 

The pressure profile is determined by the requirement 
that the base atmosphere (ie., ignoring the density fluc- 
tuations) be in hydrostatic equilibrium. The gravita- 
tional acceleration and gas density profile are sufficiently 
complex that the pressure profile is derived numerically 
and is shown in Figure [2j The core value of the pres- 
sure is Po = 4.0 X 10~^^ dyne cm~^. The pressure and 
density profiles can be combined to give a temperature 
profile (Figure |2|, which is used as mentioned to provide 
the normalization of gravity. The temperatu re profile 



was selected to resemble tvpical clusters {e.g., [Vikhlinin 
et al. 2006). As Figure [2] shows, the atmosphere has a 



core temperature Tc ~ 2.5 keV (assuming a mean molec- 
ular weight p = 0.5), with a total temperature range 
of ~ ±0.2 Tc over the computational grid. The sound 
speed, with a profile identical to that of the tempera- 
ture, features a core value cq = 895 km s~^. 

The ambient magnetic field (see Figure |3| is designed 
to be a tangled, yet analytically specifiable, structure of 
the following form: 



where the two components are 
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Fig. 4. — Average values of lOOo; = lOOPs/Pgas (solid line), 
plotted as a function of radius from the cluster center. The dotted 
line, shown for comparison, illustrates lOOo; = 1, corresponding to 
f3 ^ 100. 



Fi{r) • m _ . 
Bq — sm6'cos(m(^ 



(6) 



and 



Ba, = sm(nO) - sin(m(^) sin(2l9). (7) 

r r 

The constants m and n are chosen to ahow the field 
to vary in and 0. In our simulations m = n = 3, 
to allow for several full field reversals over either an- 
gle. Fi{r) and F2{r) are, in general, arbitrary functions 
of r that we use to approximate a constant /3 atmo- 
sphere with fiuctuations. In practice, we choose Fi{r) = 
rBa{r)[l^sin{akrir)] andF2(r) = rBa{r)[l-\-sin{akr2r)] 
to provide field variation with radius. The constants 
akri = 27r/100 kpc""*^ and akr2 = 27r/58 kpc""*^ are used 
to provide fields that vary on spatial scales comparable to 
or slightly larger than the width of the jet. These scales 
should be interpreted as the maximum length scales over 
which the field varies in our clusters rather than the 
scales on which magnetic power is concentrated. The 
strongest fields in our clusters vary over a few tens of kilo- 
pa rsecs, closer to the valu es inferred for observed clusters 
in|Vogt fc Enfilin' (2005), for example. Ba{r) is simply 
the ambient field strength as derived from the pressure 
profile and the requirement that f3 ^ 100 over a radial 
average, as shown in Figure [4j 

3. DISCUSSION 

3.1. Energy Flows 

We begin our discussion with an examination of en- 
ergy flows in these systems. By carefully tracking the jet 
power over time and comparing the integrated power to 
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Fig. 5. — Energy flow in the ST model. Upper left: a comparison 
of the known added kinetic energy (dashed line) to the measured 
change (relative to the initial value) in kinetic energy on the grid 
(dotted line) as they vary in time. The total (kinetic + thermal 
+ magnetic+gravitational) inflow energy is shown (solid line) as 
a reference. Upper right: same for the thermal energy. Lower 
left: same for the magnetic energy. Lower right: same for the 
gravitational energy. 

the measured change in energy on the grid for a given 
energy type, we can determine how energy is converted 
as the flow evolves. Additionally, we can use the passive 
color tracer Cj to indicate what fraction of this added 
energy enters the ICM. We are particularly interested 
here in the amount of thermal energy transferred to the 
ICM since this energy is most likely to be of relevance 
to the quenching of cooling flows. Again, we neglect ra- 
diative cooling; our primary goal is to understand energy 
transfer from the jets. 

Figures [5] - |8] illustrate the energy flows for each of 
the four jet simulation models we use. In each flgure, 
four plots are shown, corresponding to the four types of 
energy present: kinetic, thermal, magnetic, and gravi- 
tational. In each plot, the solid line represents the to- 
tal amount of energy added to the grid by the jets, mi- 
nus the energy measured to have exited the grid through 
the outer boundaries. In these models, the latter quan- 
tity is very small since the systems would be in hydro- 
static equilibrium were it not for the evolution of the 
added density perturbations and magnetic fleld fluctua- 
tions. In practice, the amount of energy flux across the 
outer boundaries is further reduced by the fact that the 
density, pressure, and magnetic fleld are at their weak- 
est in these regions. The dashed line in each plot shows 
the integrated jet power for that particular energy type. 
The dotted line shows the measured change in that en- 
ergy type on the computational grid, excluding the jet 
inflow region that we treat as a boundary. 

Since the systems were not in perfect equilibrium ini- 
tially, the ambient medium relaxes as the simulations 
progress, causing some exchange of energy independent 
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Fig. 6. — Same as Figure [s] but for the 126 model. 
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Fig. 8. — Same as Figure [s] but for the RE model. 
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Fig. 7. — Same as Figure [s] but for the 113 model. 
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Fig. 9. — Same as Figure [s] but for the AM model. In this case, 
no jet inflow is shown since tne jets are never activated. 



of jet action. Really, we are interested in energy flows re- 
sulting from the jet-ICM interactions rather than those 
established by the details of the ICM initialization, so 
this background relaxation is subtracted away from the 
measured change in energy {dotted lines) in each of Fig- 
ures [5]- [S] To calibrate this background energy flow, we 
used the results of the AM model, which featured no jets. 



Figure 
The so 



9^ shows the energy flow evolution for this model, 
id line in each plot in Figure [9] illustrates the ac- 



curacy with which we account for energy changes on the 
grid (the expectation value is zero) , while the dotted lines 
show the change in each energy type as usual. No dashed 
lines are included since no jet inflow was introduced. As 
expected, the AM model relaxes by gaining kinetic and 
thermal energy at the expense of gravitational and mag- 
netic energy. The amounts of energy that change form 
in the AM model are only approximately 1% of the total 
added jet energy in the ST model, but this amount is 
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Fig. 10. — Energy deposited by jets into the ICM over time, as 
a fraction of total energy added to the grid. TE: thermal, KE: 
kinetic, GE: gravitational, BE: magnetic. The model name is in- 
dicated in the upper-right portion of each plot. 

worth tracking since it exceeds the total added jet mag- 
netic energy in that model, for example. 

Focusing now on a discussion of Figures |5]-[8j we point 
out some general features of energy flow in these sim- 
ulations. Unsurprisingly, the inflowing jet energy in all 
cases is dominated by the kinetic energy component (il- 
lustrated by the dashed line beingjust below the solid 
line in the upper- left plot in Figurepj for example). Most 
of the remainder of the energy added by the jets is in the 
form of thermal energy, with relatively minor magnetic 
and gravitational components. The gravitational com- 
ponent {dashed line in lower-right plot)^ in particular, 
is roughly constant for active jets, representing the in- 
troduction of jet material in the gravitational potential 
upon jet initialization. We note that the thermal compo- 
nent of the inflowing jet power {dashed line in upper-right 
plot) asymptotically reaches a somewhat larger fraction 
of the total jet luminosity for all active jets as time goes 
on. This is caused by the sides of the jet cylinder al- 
lowing material back across that boundary, and is thus 
an artifact of the way in which we launch our jets. The 
calculated jet power is reduced by this amount, however, 
so that the depictions of jet power are an accurate rep- 
resentation of the net inflow from the jet region. 

There are also several global trends present in the mea- 
sured changes on the grid {dashed lines) for the jet-driven 
models. First, we note that in all models the measured 
thermal energy exceeds any of the other components. Re- 
gardless of the details, all models generate substantially 
more thermal energy than is introduced through the jets. 
Additionally, we note that all models feature appreciable 
increa ses in the gravitatio nal energy. Unlike our previous 
work ( O'Neill et a/.|2Q()5 ), these simulations include a re- 
alistic rnoHer"or"gmvity^ which the total gravitational 
energy on the grid is comparable to the total thermal 



energy. Since little gravitational energy is introduced 
at the jet orifice, the increase in gravitational energy is 
almost entirely a result of energy exchange in these sys- 
tems. Finally, the measured magnetic energy evolution 
also features an interesting shape, suddenly increasing in 
each figure. This is representative of the fact that the 
measured magnetic energy change on the grid is actu- 
ally negative at early times and so is not shown on these 
plots; a point we will discuss shortly. 

Before we discuss the details of each model, we should 
first introduce another measure of energy flow in these 
systems. By incorporating the passive color tracer, we 
can determine what fraction of the measured change in 
energy on the grid takes place in the ICM. This is ac- 
complished by measuring (1 — Cj)AE^ where Cj = 1 for 
jet material and Cj = for the ICM. Figure 10 shows, 
for the four jet models, this fraction for eacnaifferent 
type of energy. Again, the background change in energy 
due to the relaxation of the ICM is removed so that the 
included changes are ultimately the result of jet-ICM in- 
teractions. The most significant trend to note is that all 
models are efficient at trans ferring energy fr om jets to the 
ICM. As we found in O'Neih et aL\ (|2005|, efficient en- 
ergy transfer seems to be a general characteristic of these 
systems, and all jet models have transmitted ~ 60% or 
more of their energy to the ICM by the end of the sim- 
ulations. Significantly, thermal energy always represents 
the largest fraction of energy added to the ICM for all 
models. Unsurprisingly, given the large /3 in the jet in- 
flow, the magnetic energy always represents the smallest 
fraction of energy added. The exact distributions of en- 
ergy vary according to model, however, so we will address 
them individually. 

Focusing for the moment on steady jets, represented 
by the ST model, there are several interesting feat ures t o 
note. Similar to what we found in [O'Neill et al. ( 2005[ ), 
the active jets efficiently convert kinetic to thernial en- 
ergy. The measured thermal energy is asymptotically a 
factor of 2 larger than thermal energy added by the jet, 
which is similar to the amount of kinetic energy lost from 
the jet inflow. From the upper- left panel of Figure [Toj we 
see that ~ 40% of the total energy added to the gridends 
up as thermal energy in the ambient medium, suggesting 
that these steady jets efficiently heat their environments. 
Although it is clear that this energy is being added to the 
ICM, we also seek to discover where in the ICM this en- 
ergy appear s. As Vernaleo fc Reynolds ( 2006 ) and HeinzJ 



et al. ( 2006[ ) demonstrated, the location of deposited en-" 
ergy is at least as significant as the amount of energy 



deposited for a system to offset cooling. Figure 11 shows 
illustrations of the entropy density s = T j p^l'^ for the ST 
and 113 models. Specifically, we apply the color tracer 
and show As = s — sq for Cj < 0.10, which removes 
the initial ambient entropy sq and focuses only on en- 
hancements in the ICM. Looking at the upper-left panel, 
most of the visible As forms a sheath around the jet. 
This is caused by a mixture of jet-ICM mixing at the 
boundary (jets have low density, so at a given pressure, 
high entropy) and ICM entropy increases resulting from 
shock heating. In the saturated image {upper right), we 
can see the effects of the bow shock, although the shock 
discontinuity itself is not visible. These images suggest 
that much of the irreversibly added energy measured in 
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Fig. 11. — Color-selected entropy density {T/ p^^^) cross sections for the ST {upper panels) and 113 {lower panels) models at comparable 
jet lengths. Entropy is only shown for material that is at least 90% ICM by mass, and the initial entropy is subtracted off. Higher intensities 
correspond to higher entropies, although the intensities in the right two panels are saturated to permit visualization of the influence of the 
bow shock. 



the ICM in the ST model probably is transferred and 
remains near the jets themselves rather than uniformly 
filling the shocked ICM region. 

Figure[5]also shows that the measured gravitational en- 
ergy change represents an increasing fraction over time 
of the total added energy for steady jets, although it still 
comprises less than 10% of the total energy by the end 
of the simulation. Since the jets are light with respect 
to the ambient medium, they contain a very small frac- 
tion of the total gravitational energy. In fact, you would 
have to lift jet material 100 times higher than core ICM 
material for the same energy gain. Rather, this mea- 
sured change is caused by the jets driving ICM material 
higher into the gravitational potential. Figure |To] also 
supports this interpretation since the fraction oigravi- 
tational energy added to the ambient medium (~ 5%) 
is comparable to the total fraction estimated from Fig- 
ure |5] The increase in gravitational energy can either be 
accomplished mechanically — through lifting, dragging. 



or entrainment — or through the ICM reacting to the 
added heat and local pressure. In the ST model, drag- 
ging (ie., as in a wake) should not play a role since the 
jet is never disconnected from the source. Since the pas- 
sive variable appears to have a sharp boundary in this 
model (which you can see from the mass-fraction-selected 
entropy in Figure 11), the ICM is also not significantly 
entrained. Any lilted ICM material would have to be 
located between the jet/ICM contact discontinuity and 
the bow shock. In the direction of jet propagation for an 
active jet, this is a very small region and so not likely to 
contribute much. Since we know that the ICM is heated, 
this is the most likely source of the increase in gravita- 
tional energy for this model. 

The evolution of the magnetic energy in the ST model 
is interesting since magnetic energy is initially lost from 
the system and later increases. In fact, this is a general 
characteristic of all four jet models. Although the to- 
tal magnetic energy budget is tiny with respect to the 
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Fig. 12. — Length evolution of our model jets. Solid lines show 
the maximum extent in the x-direction as measured by the passive 
variable. Dotted lines indicate the position of the bow shock as 
a function of time. Values are averaged over the two jets in each 
simulation. 



other forms of energy (making the accurate subtraction 
of the ICM more crucial), we stih want to understand 
how magnetic fields evolve, particularly since these fields 
are an essential component of the synchrotron emission 
we observe. There are two possible physical scenarios 
one can imagine to achieve a net decrease in magnetic en- 
ergy after the ICM relaxation has been subtracted. One 
is that the stronger central ICM fields get advected up- 
ward to where they can more easily expand, reducing 
their magnetic energy. This seems implausible since the 
ST model does not advect much ambient material with- 
out compressing it. Also, the observed effect is almost 
immediate, and so operates at early times. The second 
possibility is that the ambient fields reconnect with the 
jet fields and then relax. Although reconnection as a 
process would not likely lead directly to a significant loss 
of magnetic energy, it has the potential to create new 
field geometries that could facilitate magnetic relaxation. 
This is the more plausible of the two scenarios since the 
jet fields certainly reconnect with the surrounding ICM 
(see the discussion in Section 3.3). Unfortunately, it is 
less obvious that this process would necessarily lead to 
a net reduction in magnetic energy since the stretching 
and twisting of reconnected fields could also amplify the 
fields (see, for example. Figure 18). After this initial dip, 
however, it is clear that the overall magnetic energy in 
the system is increasing. Some of this is compression 
of the ICM magnetic field by the jet/cocoon/bow shock, 
but Figure [M] suggests that much of it is simply Poynting 
fiux introduced by the jet infiow. The animations asso- 
ciated with Figure [M] further suggest that the magnetic 
energy in the cocoon is replenished by the jet since the 
field in the cocoon is reduced in magnitude during peri- 
ods of jet quiescence, presumably as a result of adiabatic 



Fig. 13. — A comparison of the magnetic energy evolution of 
the ST model to the ICM magnetic field structure. The solid line, 
corresponding to the values on the y-axis, shows the change in 
magnetic energy on the grid in the ST model. The dotted line 
(scaled for comparison) shows the magnitude of the ICM field in 
the direction of jet propagation. 

expansion. 

Interestingly, there is a correlation between changes in 
magnetic energy and jet encounters with the structure 
of the ICM magnetic field. As Figure |4] shows, the ra- 
dial distribution of magnetic field is not uniform, but 
rather peaks in magnitude every 50 kpc or so. Using the 
measured jet lengt hs i n Figure 12, which we will discuss 
further in Section 13. 21 we can translate the time evolu- 



tion of magnetic energy to evolution as a function of jet 
length. Figure 13 shows the relationship between the 
change in magnetic energy and the ICM magnetic field 
in the direction of jet propagation. The two quantities 
appear to be correlated — or anti-correlated — with a 
slight offset. The causal reasons for this are interesting 
since it is not obvious that the structure of the globally 
weak ICM field should be detectable in the energy fiow. 
Animations associated with Figure 14 show the magnetic 
field evolution for the ST and 113 models. In each case, it 
is clear that the ICM magnetic energy is enhanced each 
time the jet passes through a strong magnetic pressure 
gradient in the ICM. These waves of enhancement lag 
slightly behind the jet intersection of a high-field region, 
suggesting that the solid line in Figure [l3] is correlated 
with, but lags slightly behind, the dotted line. Since a 
substantial fraction of the magnetic energy on the grid 
is located in these regions, it makes sense that their per- 
turbations would manifest themselves in the measured 
magnetic energy changes. 

We turn now to a discussion of energetics in the two in- 
termittent jet models, 126 and 113. As Figure [12] and the 
animations associated with Figure [l] illustrate, the ap- 
pearances of the two intermittent jet models deviate from 
one another very early on. This is due mostly to the fact 
that the 126 model manages to form a well-defined low- 
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Fig. 14. — Volume rendering of the magnetic field strength for the ST (left) and 113 (right) models. Lighter colors correspond to stronger 
fields. Animations of these quantities as seen from several different angles will be available through the electronic version of this paper. 



density jet channel before the first 'off' cycle, whereas 
the 113 model must struggle at first to break through 
the dense ICM. As Figures |6] and [7| show, however, the 
two models are energetically quite similar. The kinetic 
and thermal energy influxes closely follow the level of jet 
activity. As an example, the 126 jet is on for 26 Myr, cor- 
responding to a monotonic increase in added kinetic and 
thermal energy, just like the ST model. When it switches 
off from t ~ 26 — 29 Myr, the total added energy levels 
off while kinetic energy briefly continues to convert to 
thermal at the terminal shock. While the jet remains off 
from t ~ 29 — 55 Myr, the system loses kinetic energy and 
gains gravitational energy. As the jet finishes restarting 
by t ~ 60 Myr, the cycle begins again, and similarly for 
the 113 model. It is worth noting that the gravitational 
energy increases monotonically through the entire life- 
time of the model. This is consistent both with the ICM 
reacting to the added heat (as in the ST model) and the 
dragging up of ICM material in the wakes of the unpow- 
ered, momentum-driven jet plumes. This latter process 
almost certainly takes place in these models since we see 
the kinetic energy of the system decrease while the grav- 
itational energy increases at t ~ 29 — 55 Myr. The frac- 
tional energy plots (Figure 10) show that both models 
are very efficient at heating tHe ICM, although the frac- 
tion of energy in the ICM at a given time can vary by as 
much as ~ 20%, depending on the level of jet activity. 

One interesting point concerning the 113 model is that 
there are enough cycles of jet activity that we can at- 
tempt to discuss asymptotic behaviors. Compared to the 
ST model, the 113 model has approximately the same to- 
tal efficiency (~ 60%) in transferring energy to the ICM. 
The fraction of kinetic and thermal energy entering the 
ICM is similar for the two runs, but the gravitational 
energy fraction is ~ 10% higher for the 113 model. This 
is slightly misleading, however, since the total amount 
of gravitational energy introduced into the system at a 
given time is higher for the ST than the 113 model. The 



fractional gravitational energy differences in Figure 10 
are thus more indicative of the difference in the total en- 
ergy available to each model. Additionally, we can again 
examine the entropy to see where the intermittent jets 
characteristically deposit their energy. The bottom two 
panels in Figure [TT] show the color-selected As for the 
113 model. Although the added ICM entropy still forms 
a sheath around the jet cocoon, there is now a much more 
substantial amount of material nearer to the cluster core. 



Likewise, as we will see in Section 3.2, the intermittent 
cocoons are characteristically fatter than those of the ST 
model. This means that more high-entropy material is 
spread through a larger volume and nearer to the cluster 
core, both of which are needed to offset core cooling. 

We turn lastly to the RE model, which evolves iden- 
tically to the ST model until the jet activity is shut off 
completely at ~ 26 Myr. In Figure |8j we see the ex- 
pected inflowing and total energy plateaus appearing at 
~ 26 Myr. As in the 126 and 113 models, the system be- 
gins to exchange kinetic energy for gravitational energy 
as the thermal energy is approximately maintained. Af- 
ter t ^ 60 Myr, however, we see the total thermal energy 
start to decrease as the kinetic energy mildly increases. 
This feature is also marginally present in the 126 model, 
but in that case the jet reactivates before it becomes 
very pronounced. In the RE model, the slope of the 
added gravitational energy also begins to flatten. This 
behavior is consistent with the transition from a system 
actively driven by jets to one in which a passive plume 
consisting of a mix of ambient and jet /cocoon material 
is gradually slowed by the drag force associated with its 
motion through the ICM (see Section |3.2| for details). 
As the relic plume rises, it should maintam approximate 
pressure equilibrium with the surrounding ICM, expand- 
ing and losing thermal energy in the process. Likewise, 
reversible heat added to the ICM can be transformed 
back to kinetic energy through adiabatic expansion. Ul- 
timately, the total gravitational energy should decrease 
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when the system becomes dominated by buoyancy, but 
there is still sufficient momentum in the plume through 
the duration of our simulations that we have not yet 
reached that regime. Figure [lO] shows that approximately 
~ 75% of the jet energy has entered the ICM by the end 
of the RE simulation. Despite being more efficient at 
delivering what energy it has, however. Figure [S] shows 
that the RE model introduces less total energy into the 
ICM than any of the intermittent or steady models as a 
result of its lower integrated luminosity. 

3.2. Dynamics and Morphology 

Here, we present a discussion of the dynamics and mor- 
phology of our simulated flows. As mentioned in the 
previous section. Figure [l] and the associated animations 
clearly show the evolution of flow speed in each jet model 
over time from various angles. Likewise, Figure 16 and 
its associated animations show volume renderings of jet 
entropy surrounded by more transparent isosurface ren- 
derings of ICM entropy. Animations such as these are 
typically the best way to quickly visualize causal relation- 
ships between flow structures and to convey the general 
character of the flow. 

First, we discuss the general properties of jet length 
evolution as a function of time, as illustrated by Figure 
12 In that figure, the solid lines represent the maximum 
distance between the cluster core and Cj > 0.90 in the di- 
rection of jet /relic motion and averaged over the two jets 
in each model. The dotted lines show the position of the 
forwardmost shock, also in the direction of jet/relic mo- 
tion and averaged over both jets. To locate the shocks, 
we used an identifica tion strategy similar to that em- 
ployed in Ryu et al. (2003) in which the temperature 
and entropy gradients were of the same sign, the veloc- 
ity divergence was negative, and the temperature jump 
was that of a Mach 1.3 shock or stronger. In Figure [T2| 
we see that substantial separation develops between the 
bow shock and the jet material at some point in every 
non-steady jet model. When the dotted line and solid 
line are exactly coincident, this corresponds to a shock 
with M < 1.3 or no shock at all. 

There are two interesting evolutionary features in Fig- 
ure [12] First, we note that the ST, 126, and 113 mod- 
els allasymptote to similar slopes, suggesting a common 
power-law evolution oi I (x for these models, where 
(5 1.3. A simple expectation for the adva nce speed, 
gotten fro m self-similar scaling arguments in [O'Neill et| 
al. (2005), for example, is that 



(8) 



assuming a homologous structure for the cocoon and that 
the ambient density follows a power law pa oc x~'^. Al- 
though the ambient density in this case is not a strict 
power law, we can estimate z^: ~ 1.3 around r ~ 100 kpc 
and n ^ 1.7 around r ~ 300 kpc. These values produce 
an expected (5 ~ 0.8 — 0.9, which is substantially flatter 
than the measured value of ^ ~ 1.3. Self-similarity is thus 
probably not a reasonable assumption for an atmosphere 
whose density dependence changes slope so dramatically. 
That the measured values of 5 for the 126 and 113 models 
are slightly flatter than the ST model is not surprising, 
however, since the jet momentum is probably distributed 
over a slightly larger area, as we will see momentarily. 



In contrast to systems in which jets actively drive the 
flow, the advance speed of the RE model asymptotes to 
a power law of approximately 5 ~ 0.4 as it decelerates 
in the ICM. Although the residual jet flows terminate 
quite rapidly after being shut off, we expect the overall 
momentum of the system to be distributed in a mix of co- 
coon material and ICM that has been entrained, dragged, 
or otherwise pushed into motion. We can estimate the 
rate at which the plumes are decelerated, assuming that 
the dynamics are dominated by the effects of drag re- 
sulting from motion of the relics through the ICM. The 
drag force is given by Fd ~ —CdApav'^^ where C^/ ~ 1 is 
the drag coefficient, A is the cross-sectional area of the 
plume, and v is the relic velocity. We assume, as before, 
that the ambient density follows a power law, this time 
specifying it more fully as pa = Po,off(^/^off)~'^, where 
the "off" subscript refers to the value of that quantity 
at the position of the jet head when the jets are shut 
off. If we assume that the radial expansion is adiabatic, 
then PV^ (X Pr^^ is conserved, assuming the volume 
of the relic is proportional to the cube of its radius. If 
the ambient pressure profile is proportional to that of the 
density {i.e., the ICM is approximately isothermal), then 
, meaning that the relic area can be approxi- 
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mated as A ^ 7rr^ = 7rrl^{x / Xof^)'^'^^^^ . Under these 
assumptions, the relic deceleration resulting from drag 
can be written as 



(9) 



where, writing the total mass of the relic as ttIyq and 
assuming = 1, 



constant. 



(10) 



These expressions can be simplified considerably when 
we realize that, measured at x^^ > 150 kpc, n ^ 1.6 — 
1.7 7. This reduces the drag deceleration expression to 
X ^ ^x~^x^^ which admits a solution x cx To 
estimate the value of ^, we can model the initial shape of 
the relic plume as a cylinder of mass = Pre(7rr^^Xoff). 
Again using the fact that />: ~ 7, this allows us to write 
^ ^ — Pa,off/Pre- Unfortunately, using the passive color 
tracer, Cj, to estimate the relic density would fail to ac- 
count for the any contribution to the total mass from 
the ICM component of the relic. To more accurately 
estimate the relic density, we instead measure from our 
simulated data the average density of material moving 
upward with velocities above a certain threshold and po- 
sitions between the jet origin and the top of the relic 
(thus excluding the bow shock). Choosing only material 
with Vx > 1.5co, which would correspond everywhere to 
initially supersonic motion in the unperturbed ICM, we 
measure an average relic density pre ^ 0-15 at t ~ 26 
Myr, just after the jets are shut off. At that relic height 
(^off ^ 155 kpc), Pa,off ^ 0.20, resulting in an expected 
value of ^ ~ —1.3. Since / cx the expected 

value of (5 is 1/(1 — ^) 0.43, which is similar to the 
observed ^ ~ 0.4. The closeness of the predicted and ob- 
served power laws indicates that this deceleration model 
is plausible, although one also assumes that this scal- 
ing would eventually break down when the relic has lost 
enough momentum that buoyancy becomes the dominant 
force in the system. 
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Fig. 15. — Shock structures in the 113 model at early {upper left), moderate (upper right), and late (bottom left) times, illustrated by 
volume renderings of —V • v. The lower-right image shows the ST model at a late time for comparison. An animation of this quantity as 
seen from several different angles will be available through the electronic version of this paper 



Moving on to a discussion of jet and flow morphol- 
ogy, we flrst describe steady jets as exemplifled by the 
ST model. Our steady jets are narrow, stable struc- 
tures, a common feature in such simulations. Since we 
do not perturb them at the source, but instead allow the 
ICM substructure to seed instabilities, they do not de- 
ve lop complicated shock stru ctures like those obs erved 



[Tregillis et al] ([2001b|) and lO^Neill et al\ (|2QQ5|, and 
the jets are not broken apart, as m Heinz et al] ( |2QQ6 ). 
Furthermore, the jet magnetic fleld is sufficient ly weak 
that the jet does not become unstable to pinch insta- 
bilities, nor does it develop an obviously modifled nose 
cone as has been seen in some str ongly magnetized jet 
simulations {e.g., 



Clarke et al. 



1986^ . 



As shown in the 



lower-right panel ot .b igure |15[ the bow shock is a rela- 



tively smooth surface and the internal shock structures 
are regularly spaced. From the animations and image 
of the flow speed associated with Figure [ll we see that 
the jets are not completely axisymmetricTnowever, and 
these slight asymmetries become much more pronounced 
when the AGNs are shut down. The cocoons of shocked 
jet material in steady jets are relatively thin and do not 
develop a strong backflow. The cocoon structure is seen 
clearly in a volume rendering of the entropy shown in 
Figure [16) where the low-density plasma is very bright. 
The ICM pressure gradient in these models is suflicient 
to prevent much material from traveling very far opposite 
the direction of jet propagation. 

The 126 and 113 intermittent model jets feature less 
axisymmetry as the flows evolve. Figure [l] shows that 
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^^^^^^^^^^^ 




Fig. 16. — Volume rendering of entropy for the ST (upper left), 126 (upper right), 113 (lower left), and RE (lower right), after the jet 
disturbances have nearly reached the computational grid boundary. Higher entropy is shown in lighter colors, and the ICM is only shown 
in isosurfaces. Animations of these quantities as seen from several different angles will be available through the electronic version of this 
paper. 



portions of the flow with moderate speeds at late times 
tend to bend in response to the ICM substructure. Like- 
wise, we can see filaments of low-speed material that form 
in the cocoon. The shock structures in these two models 
are very different from the steady case. Figure [15] shows 
a volume rendering of the 113 shock structures at three 
different times. Not only do the 113 shocks look differ- 
ent from those in the ST model, but they look different 
at a given time from earlier or later shock structures in 
the same model. The animation associated with Fig- 
ure [15] clearly shows the propagation of internal shocks 
and waves associated with restarted jet activity. Typi- 
cally, the subsequent jets form a bow shock in the old jet 
material that quickly advances forward through the low- 



density channel. This was seen e ven in the earliest mod - 
els of restarting jets conducted by Clarke & Burns ( 1991 ). 
A set of irregularly shaped reverse shocks that dampen 
into waves moves back through this material each time 
the jet turns off. 

Compared to the ST model, the two intermittent mod- 
els feature much wider cocoons that are 'puffed out' by 
subsequent episodes of jet activity. Figure [16] shows im- 
ages of the entropy for these models. In the 126 model, 
we can see that the second episode of jet activity has 
run into the plume from the first episode, causing it 
to widen. Likewise, the 113 model cocoon has several 
visible width enhancements associated with subsequent 
jet activity. These features are associated in time with 
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Fig. 17. — Example of magnetic reconnection in the ST model. The light-colored, straight streaklines pointing from the bottom left to 
the upper right of each image represent the jet velocity, which is surrounded by a toroidal magnetic field. A set of ICM magnetic field lines 
sits in the path of the jet, color coded to magnetic field strength (light: strong). From t ~ 20 Myr (upper left) to t ~ 40 Myr (lower right), 
the jet field pushes through and eventually connects with the ICM field structure. 



the measured grid increases in thermal energy, represent- 
ing energy deposition in the plumes. Interestingly, these 
structures persist after they are formed, presumably be- 
cause subsequent jets are easily redirected by the ICM 
into the relatively low-density channels formed by previ- 
ous jets. The animations of entropy for the 113 model in 
particular show these features lasting to the end of the 
simulation. There are at least two reasons, both having 
to do with where shocks form, that intermittent jets pro- 
duce wider cocoons than their steady counterparts. The 
first is that the ICM is drawn in behind the wake of the 
rising plume before the next intermittent jet launches. 
This refilling is incomplete, as seen in the animations as- 
sociated with Figure[^ but is sufficient to cause the next 
episode of jet activity to produce a shock at the base of 
the plume. The second reason is that the restarted jet 
forms a terminal shock and multiple reverse shocks of its 
own (as seen in Figure 15) as it propagates through the 
existing plume. Both oTthese episodes of shock forma- 
tion deposit thermal energy in the plume and cause the 
cocoon to expand laterally as a result. 

Finally, the relic model is distinct from all other mod- 
els in that the jet plumes completely detach from the 
source. Figures [l] and [16] show the evolution of this 



model, in which two plumes of material rise in the ICM 
as the deposited momentum is depleted. When the jets 
shut down, the fiows are already well over 100 kpc from 
the cluster core, and their cocoons are prolate in shape. 
As the ends of the plumes decelerate, however, the relic 
bubbles grow less prolate. Since observed relic and ghost 
bubbles are typically seen to be intact structures, we note 
that the entropy at late times shows relatively rounded 
morphology, particularly near the leading edge of the 
bubble. Significantly, we do not observe the formation of 
a torus, in contrast with simulations of bubbles initial- 
ized as spherical cavities in cluster s that would o therwise 
be in hydrostatic equilibrium (see O'Neill et a/.|[2009 for 
an extensive discussion). Finally, it is important to note 
that we have only explored a small part of the avail- 
able parameter space with this particular relic model. 
To fully evaluate whether jet/ICM magnetic fields st abi- 
lize jet-blow n relic bubbles, as has been suggested by |De| 
Young' ^2003 ) , for example, further variation of magnetic 
held strength and geometry would be required. 

3.3. Magnetic Field Evolution 

Finally, we devote some space to a discussion of m ag- 
netic field structures in these fiows. Section 13.11 ad- 
dressed the evolution of magnetic energy in these sys- 
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terns, but the details of field structure evolution were not 
extensively disc ussed. Employi ng a visualization package 
called FieldVis ( [Field et a/;||2QQ7| ), we are able to graph- 
ically represent magnetic neld lines and to map useful 
flow properties on them as needed. This is a tremen- 
dously powerful approach since the challenges of visual- 
izing multiple three-dimensional vector and scalar fields 
can otherwise be overwhelming. 

Figure [Tt] shows the evolution of magnetic field struc- 
tures in the ST model. At early times {upper-left image) ^ 
the jet is propagating toward an undisturbed ICM field 
structure, similar to one of the structures seen in Fig- 
ure [3j In that image, the jet moves from the lower left 
toward the upper right of the image and is represented 
by a set of lightly colored velocity streaklines. The jet is 
surrounded by a toroidal field, represented by the darker 
set of coiled lines. As the jet disturbance first reaches the 
ICM {upper-right image) ^ the ICM field lines are pushed 
aside and stretched. The intensity of the magnetic field 
lines represents the field strength (brighter=stronger) , 
and we see that the field is enhanced as it wraps around 
the jet. This process continues in the next frame {lower 
left)^ where the jet has completely penetrated the ICM 
field. In that same image, we also note evidence of an 
internal shock within the jet corresponding to a concen- 
tration of field lines roughly halfway down the jet struc- 
ture. In the last image {lower right), we see that the jet 
field has interacted with the ambient field. To explain 
this image completely, we must first point out that the 
identified ICM field lines were selected to lie ahead of 
the jet so that the jet field would necessarily intersect 
them. This final image shows that at least one of the 
ICM field lines has become attached to the toroidal jet 
field. This is seen most easily by noting that there are 
four visible lines to the left of the jet, but only three 
emerge to the right while the other is attached to the 
larger toroidal structure near the end of the jet. We ex- 
amined the values of the passive variable in this region 
to find that in fact the large toroidal field is associated 
with jet material and connects with a line that is asso- 
ciated with the ICM. Although reconnection in our nu- 
merical scheme does not represent a significant amount 
of direct energy exchange, the geometries that develop 
from it have the potential to affect the system energet- 
ics. As this system evolves, the reconnected field line will 
get stretched out with the jet flow, possibly leading to 
an alignment between the jet and field orientation. This 
radial realignment of ICM fields could have important 
implications for cluster heating via thermal conduction 
along field lines, as has been sugg ested by [Bogdanovic 



et al. 



tiona. 
from 



(|2009j) and Quo fc Oh| ( |2QQ9| , for example. Addi- 



ly, this alignment could be similar to that inferred 

polarization maps of some stronger sources (Bridle 

et a/.||199 4). We list this outcome only as a specula- 
tive possibility, however. Constructing synthetic radio 
observations is the only rigorous way to determine the 
polarization properties of a simulated structure, and we 
defer that discussion to a later work. 

Interestingly, the ICM magnetic field seems to inter- 
act strongly with the behavior of the intermittent jets. 
Although the nominal mean magnetic pressure is a fac- 
tor = 1/100 smaller than the gas pressure, local 
magnetic pressure gradients can still exert a substantial 
force on the flow. There are two reasons for this. First, 




Fig. 18. — Sample magnetic field evolution in the RE model. 
Upper panel: the initial ICM magnetic field configuration at 
r ~ 200 kpc from the cluster core. Intensity refiects field strength 
(light=strong). Lower panel: the same region at t ^ 120 Myr. 
Now, the relic velocity is shown as streaklines with 'baseball dia- 
mond' glyphs pointing in the direction of motion. The surrounding 
stretched lines represent ICM field lines that have gotten stretched 
and amplified by the motion of the relic. 



the variations in field strength can lead to locally much 
smaller values of (3. Second, it is important to remem- 
ber that the dynamical consequence of magnetic pres- 
sure comes through its gradient, not its magnitude, per 
se. When the local length scale of magnetic field varia- 
tions is smaller than the length scale of the gas pressure 
variations, the dynamical role of the magnetic pressure is 
stronger than the (3 parameter implies. In the ST model, 
these forces are swamped by the momentum of the jets 
and are not easily seen. In the intermittent cases, how- 
ever, there is some evidence that the ICM field influences 
the cocoon shape as the jets propagate. Figure [M] and 
the associated animation illustrate this point well.T^ the 
jet restarts, the bright regions of enhanced field in the jet 
coincide with the ICM enhancements. Likewise, we can 
see that the cocoon width responds to these fields. Since 
this is more difficult to see in other models, however, we 
conclude that it is the interaction between intermittency 
and the ambient field that makes these structures stand 
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out. 

Figure 18 shows another field geometry, this time for 
the RE model. The upper image shows a set of undis- 
turbed ICM field lines at roughly 200 kpc from the clus- 
ter core. These have a structure much like those seen in 
Figure [3] and the upper-left panel of Figure 17 Again, 
field strength is shown as a surface property of the lines 
and the structure of the constant-/^ atmosphere can be 
seen from the darker background fields. The lower image 
shows the interaction of the relic bubble and the ICM 
field, where the relic is represented by velocity streak- 
lines using a different glyph. As the relic moves from the 
lower left to upper right of the image, it drags a signif- 
icant amount of magnetic field with it. Specifically, we 
see several very strong (i.e., brightly colored) field struc- 
tures that are stretched in the direction of relic propaga- 
tion. These are ICM fields that connected with the relic 
bubble fields as it was rising. Once this happens, the 
continued propagation of the relic serves to stretch and 
amplify these fields. Again, the passive variable confirms 
this explanation by showing that several of the strong 
field lines in the lower image are associated with both 
jet and ICM material. This is a prime example of how 
magnetic fields can get enhanced and transported in a 
cluster even when the central AGN engine is off. 

4. CONCLUSIONS AND ASTROPHYSICAL IMPLICATIONS 

We have conducted a series of simulations designed to 
examine the role of intermittency and cluster substruc- 
ture on the energetics, dynamics, and morphologies of 
jet- driven fiows. 

The most important results from this work are sum- 
marized here as follows. 

1. All models of jet behavior, including steady, inter- 
mittent, and relic flows, are at least ~ 60% eflicient at 
delivering thermal energy to their environments. While 
this efliciency is remarkable, none of the jet models are 
particularly adept at distributing this energy uniformly 
throughout their environments. Specifically, the transfer 
of energy from steady jets to the ICM takes place very 
near the jet cocoon, while intermittent jets spread the 
energy in only a slightly wider pattern. Jet-blown relic 
bubbles are the most eflicient at transferring what energy 
they have, but they also transmit less total energy than 
any other model. Taken together, these results confirm 
that the lack of u niform heating in h ydrodynamic sys- 
tems described by Vernaleo & Reynolds (2006) is also a 
problem for MHD Jets in realistic magnetized clusters. It 
is possible t hat the subgrid modeling of tur bulence un- 
dertaken by 
to alleviate t 



Briiggen fc Scannapieco (2009) could help 
lis problem in magnetized systems, too, but 
this will require a modification of the subgrid model that 
has been used thus far. 

2. The morphologies of intermittent jets retain evi- 
dence of previous episodes of jet activity. Specifically, 
subsequent bouts of jet activity encounter previously 
formed structures, creating reverse shocks at their in- 
terfaces. The cocoons formed by such fiows also evolve 
to be characteristically wider than those of steady jets. 
These structures can be associated with particular depo- 
sitions of energy into the ICM and may indicate where 
heating is taking place. A full exploration of whether 
such features are observationally detectable is deferred 
to the analysis of synthetic observations from our simu- 



lated data ( Mendygral et a/]|in prep ). 

3. The presence of cluster substructure can perturb ini- 
tially axisymmetric jets and does so much more quickly 
when the jets are non-steady. This is a very natural 
way to perturb jets in clusters that have not undergone 
recent mergers and which appear to be relaxed. Still, 
as noted, the density and magnetic structures included 
in these simulations are small and insuflicient to signifi- 
cantly isotropize the energy distribution in cluster cores. 

4. Magnetic structures in the ICM are affected by re- 
connection and can get carried along with jet-driven and 
passively evolving fiows. Although the specific process 
of reconnection will never directly increase the magnetic 
energy, it can facilitate field geometries that operate in 
conjunction with the fiow to amplify fields through field 
line stretching and twisting. Additionally, such processes 
have the potential to advect core ICM fields into the 
outer regions of clusters and could possibly orient such 
fields in the direction of jet motion, with potentially im- 
portant implications for heat transfer via thermal con- 
duction. Interestingly, local ICM magnetic pressure gra- 
dients are capable of affecting these fiows even when the 
global field strengths are considered to be insignificant. 
While the magnetic energy increases associated with such 
interactions are minor, this is a specific way in which 
MHD jet fiows with nominally weak fields differ from 
their hydrodynamic counterparts. 

5. The motion of relic plumes produced by AGN jets is 
initially determined by a combination of the momentum 
resident in these systems and the drag force induced by 
moving through the ICM. In our model, this mode of evo- 
lution persisted for 100 Myr after the jet had switched 
off, suggesting that buoyant models of jet propagation 
in clusters may only be relevant at late times or for jets 
with less total momentum deposition. The resulting relic 
morphology and evolution is a useful reminder of the 
important differences between simulations of jet-formed 
bubbles and bubbles that are initialized as spherical cav- 
ities. 
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APPENDIX 

APPENDIX: ESTIMATING THE JET LUMINOSITY 



As noted in Section |2.1[ the ratio of the unit kinetic to thermal luminosities contains a prefactor that is of order 
unity. To see how this prefactor depends upon the details of jet injection, let us first define our jet density, pressure, 
and velocity profiles: 



Pj = Pi,o 
Pj = Pj,o 



r <rs 
r < Ts 

r < Tn 



(Al) 
(A2) 

(A3) 



where Vj is the outer radius of the jet core, (> Vj) is the outer radius of the transition sheath, and the "zero" 
subscript indicates a baseline value that could in principle be time-varying. To understand how these profiles affect 
the luminosity, we first compute the kinetic luminosity, L/^: 



1 







PjV^dA- 



^Pj,OVj,0 



rdr 



(A4) 



20 "^'^ ^ 

Following a similar procedure for the thermal luminosity, L^: 

Lt 



PjfiVjfiif^ + 3/i + 6); 11 = Ts/rj. 
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\ ^ p] 


Jo 


[(7-1) d 



VjdA 



7rr|7 
3(7-1) 



(A5) 



Pj,oVj,o{p'^ + 1) 



We see that each of the two luminosities has a dependence upon the size of the transition region as parameterized 
by /i. That the two luminosities vary differently with ja simply refiects the decision to have the velocity transition 
smoothly to zero, independent of the density or pressure. To recover the prefactor found in Section |2.1[ we take the 
ratio of the two luminosities: 

Lk ^ 3(7-l)p,,o^|,o(/i'+3/i + 6) 

Lt 207p,-o(/i2 + /i + l) ^ ^ 

Plugging in the values from Table 1, pj^o = 8.33 x 10~^^ g cm~^, pj^o = 4.00 x 10~^^ dyne cm~^, and Vj^o = 2.69 x 10^ 
cm s~^, we obtain: 



L~t 



0.904 



P 



3/i + 6 



/i + 1 



(AT) 



For our assumed jet profile, p = Vs/vj = 7 kpc/3 kpc, so the ratio of kinetic to thermal luminosity is L^/Lt = 1.90. 
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